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We investigate the origin of the long period oscillation of 
the interlayer exchange coupling in Fe/Cr trilayer systems. 
Within the stationary phase approximation the periods of the 
oscillations are associated with extremal vectors of the Fermi 
sphere of Cr. Using a realistic tight-binding model with spin- 
orbit interaction we calculate the coupling strength for each 
extremal vector based on the spin-asymmetry of the reflection 
amplitude for a propagating state impinging from the Cr to 
Fe layer. We find that for (001) and (110) growth direction 
the biggest coupling strength comes from the extremal vector 
centered at the ellipsoid N of the Fermi surface of Cr. 



I. INTRODUCTION 

The magnetic properties of multilayer systems have 
been a subject of intensive study, particularly for the 
last decade. The initial observation that two ferromag- 
netic layers can be coupled antiferromagnetically, when 
separated by certain type of magnetic or non-magnetic 
spacer (Griinberg et alu), and the following discovery 
that this coupling possesses a damping oscillatory be- 
havior. (Parkin et al.a) stimulated a great deal of in- 
terests. Moreover, because, o£ the phenomenon of giant 
magnetoresistance (GMR)q'EI, these structures have be- 
come very promising candidates for applications in the 
magnetic recording industry. 

The construction of ultrathin structures poses great 
challenges from the experimental point of view. The ini- 
tial samples grown with sputtering displayed the char- 
acteristic oscillatory behavior, however, soon after the 
first measurements, the demand for strict control over 
stoichiometry and disorder made the growth of excellent 
epitaxial sandwiches or superlattices necessary. In the 
case of a sandwich, it is a common practice how that 
one uses a wedge geometry for the spacer layeiu. In this 
way one can study the spacer thickness dependence of 
the interlayer exchange coupling (IEC) using only one 
sample. 

From the theoretical perspective, the calculation of 
IEC has also attracted a lot of attention. Several differ- 
ent approaches to the problem have been suggested and 
used extensively. One can divide all these approaches 
into two major classes. One class consists of total energy 
calculations, in which one calculates the difference in the 
total energy of the sample for the two important magne- 
tization configurations (the magnetization of the ferro- 
magnetic layers being parallel in one case, antiparallel in 



the other). Local spin density approximation (LSDA) 



(see for example R-efs.. 
binding (TB) methodsot 



and semiempirical tight- 
have been used for this pur- 
pose. The second class consists of various model calcu- 
lations. In this we include the calculations based-atuthc 
perturbative treatment of the RKKY intep|Ctionij'E£l, an 
adapted Anderson (or sd-mixing) model! 1 ii 14 !, and per- 
haps most importantly, we include all the various cal- 
culations of IEC that exploit the idea of quantum in- 
terference and quantum confinement in the spacer ma- 
teriall!5li3. The multiple reflections that the electrons 
experience at the spacer/magnetic layer interfaces bear 
an analogy to Fabry-Perot like interferometry. Our own 
method belongs to this class. All the model calculations 
have certain approximations built in, and hence, one has 
to be aware that in some cases the accuracy of the re- 
sults is limited. Nevertheless, the potential for numer- 
ical errors involved in these calculations is significantly 
smaller than in the total energy calculations, particu- 
larly for large spacer thickness, and the computational 
effort is much less as well. 

It is now widely accepted that the periods of the oscil- 
lations of IEC are given by the extremal vectors of the 
Fermi surface of the spacer material. In this paper we 
are going to study IEC of Fe/Cr/Fe, so the spacer is Cr. 
Chromium has a complicated Fermi surface with a rich 
variety of spanning features. Therefore one would ex- 
pect in principle a multiperiodic oscillation of IEC with 
respect to the Cr layer thickness. However, only two pe- 
riods have been observed in the experiments so far. The 
short period is roughly equal to two monolayers (ML) 
of chromium and it is believed to be due to the nesting 
vector that gives rise to antiferromagnetism in Cr. The 
origin of the long period (measured to be somewhere be- 
tween 15 A and 18 A) has been proven a more difficult 
problem. 

We use an empirical TB method in order to calculate 
the coupling strengths associated with different periods. 
By using the "force theorem"!!!!, one can find the change 
in the density of states for the trilayer system with re- 
spect to the bulk case. Electrons with different spin see 
different effective potentials at the Fe/Cr interface, and 
consequently, the multiple reflections are sensitive to the 
magnetizations of the Fe layers being parallel or antipar- 
allel. The base of the analytical work was developed irua 
series of papers by Brunoli3. Stilesta and SlonczewskiEa 
also showed how one can calculate IEC based on the 
spin asymmetry of the reflection amplitudes. In more 
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recent papers the method was applied with encouraging 
results in real systemsEj. We have used the same method 
as Lee and Chang used for Co/Cu/Co systemsEil, but 
with some modifications, in order to take care of the 
much more complicated Fermi surface of Cr. The main 
conclusion of our work is that the long period oscilla- 
tion originates from the extremal vector at the ellipsoid 
centered at point N of the Cr Fermi surface. This is 
true for both (001) and (110) growth directions. The 
result agjeees with that obtained with LDA calculations 
by StilestEl. We will also discuss its validity in v,iew of 
the outcome of recent photoemission experimentscj. 



II. THEORETICAL MODEL 



The Interlayer Exchange Coupling is defined as 
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where f2 f and CIaf are the grand canonical potentials for 
the ferromagnetic (F) and antiferromagnetic (AF) con- 
figurations and S is the area of the sample. The electrons 
are partially reflected at the interfaces of the Fe/Cr/Fe 
sandwich, and the change in the grand canonical poten- 
tial is given by the force theorem (or the frozen potential 
approximation) as 

i r + °° 
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where v labels the two configurations (F and AF), Go is 
the bulk Green's function of the spacer material, T is the 
T-matrix, / is the Fermi-Dirac distribution , and L and 
R stand for the left and right interface respectively. Tr 
denotes the trace over the z-component of the wavevector 
and the spin index. The asymptotic form of IEC (for 
large thickness D of the Cr layer) is given in first order 
by the simple analytical result (see Ref. |2l] for details) 
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In Eq. (3), qfj — k" — kf is an extremal vector of 
the Fermi surface of Cr parallel to the growth direction 
(which is taken to be the z direction). The superscript 
a labels different kjj 1 , that is the in-plane component 
of the wavevector. kf and kf, are the z-components of 
the wavevector of the incident (reflected) and reflected 
(incident) electron for the right (left) interface, and uf 
is the group velocity at a point (k^,kf) of the Fermi 
sphere. The index a labels the two spin mixed states 
that are degenerate at each (kfpfcf) (Kramers degener- 
acy). The reflection amplitude from a state (kf ,a) to a 
state (fc", a') is denoted by rf^ v (where I = L,R labels 
the interfaces), nfj is related to the curvature radii at 
the two endpoints of qfj, 
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is called the geometrical weight of the extremal vector 
qfj. D is the thickness of the Cr layer, #y is a con- 
straint function, which is 1 for extremal vectors and 
otherwise. The phase factor <pfj ~ 0, tt/2, and it for 
maximum, saddle, and minimum extremal points respec- 
tively. We perform the calculation at zero temperature, 
so F£(D,T) = 1. For finite temperatures F%(D,T) is 
less than 1, if the spacer is magnetidl3. 

The reflection amplitude factor Ai?^ depends explic- 
itly on the reflection amplitudes, which, as we mentioned 
above, are calculated from an empirical TB model with 
s, p 3 , and d 5 orbitals. Because of the sandwich geome- 
try one needs to include, not only the propagating Bloch 
solutions, but also the evanescent states, i.e. the states 
with complex wavevectors. The details of the method 
can be found in Ref. ^4]. For Cr, the spin-orbit coupling 
is significant only for the so called lens area and at the 
point across the TH line where the electron surface al- 
most touches the hole surface (see Fig. [l]). For all other 
vectors the spin-orbit interaction does not play an im- 
portant role, the spin is approximately a good quantum 
number, and one finds that 
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where the superscript + (— ) denotes majority (minority) 
(5) spin electron. 



and 
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III. Fe/Cr/Fe SANDWICHES 

As we pointed out in Sec. I, the Fe/Cr/Fe case has 
been a difficult one for the calculation of IEC, and this 
is mainly because of the complicated Fermi surface of Cr. 
For our own TB calculations we have used the parame- 
ters from Ref. These parameters do not include any 
relativistic effects, and particularly spin-orbit coupling. 
We have included spin-orbit coupling with its character- 
istic parameter £ = 2.84 mRy taken from Ref. ^6|. This 
value agrees well also with the value found in Ref. |l7| . 
The (100) cross section of the Fermi surface, that we 
obtained, is shown in Fig. [|. In this figure we used a 
slightly larger parameter £, so that the effects of the 
spin orbit interaction become easier to see. The Fermi 
surface we obtained (with the correct £) agrees very well 
qualitatively, and in most cases quantitatively, with the 
ones obtained from experiments and other more involved 
band structure calculations (see for example Ref. |27| ). 
The inclusion of spin-orbit coupling affects significantly 
the area around the lens vector, as we mentioned before. 
One needs to take this effect into consideration, because 
the period related to this vector is very close to the long 
period of the oscillations of IEC. 

Our results are summarized in Tables || and [n] for the 
two growth directions that we have studied. In the ta- 
bles, c = |, ^= is the distance between two Cr layers for 

(001) and (110) orientation respectively (a = 2.88 A is 
the Cr lattice constant), and A 's are the periods that we 
found. In the last column we give the coupling strength 
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associated with each vector q°j. Gfj is the geometrical 
weight and n is the number of equivalent extremal vec- 
tors in the first Brillouin Zone. For the thickness D we 
take a typical value of 10 A. One can see that |Ai?| and 
vary significantly for different extremal vectors, and 
therefore, the overall coupling strength is not negligible 
only for very few vectors. In the next two subsections we 
shall discuss the analyses for (001) and (110) orientations 
separately. 



A. (001) orientation 

In Fig. [I], we give the (100) cross section of the Cr 
Fermi sphere. The wavevectors are measured in — . We 
have marked some extremal vectors, which, according to 
the theory, are candidates for the observed periods of the 
oscillations of IEC. Four vectors are of particular impor- 
tance. One is the vector (Ni) at the ellipsoid pocket at 
N. Two other special vectors are the one spanning the 
lens (Li) and the one just outside the lens (£2)- Fi- 
nally, we have the nesting vector (V n ), which connects 



two almost parallel lines, one from the electron octahe- 
dron around point T, and the other from the hole octa- 
hedron around point H. 

For the extremal vector N±, both Gfj and |Ai?| are big, 
and hence, the coupling strength is large (7.94 mJ/m 2 ). 
The next closest amplitude is 2.80 mJ/m 2 , for the ex- 
tremal vector X that spans the electron ball at point 
X [k|| = (0,0.42)^]. This is still 2.83 times smaller 
than the ellipsoid contribution. All the other ampli- 
tudes, with the exception of the nesting vector, are, ap- 
proximately, at least one order of magnitude smaller. In 
Figs. ||(a)-(b), we present the behavior of the modulus 
of the reflection amplitude \r + \ (|r~|) of the majority 
(minority) spin electron as a function of the energy for 
k,| = (0,0.5)^ and k y = (0,0.42)^. For these ley's, 
the spin-orbit coupling is not very important and |Ai?| 
is given explicitly by the square of the difference of r + 
and r~. These two complex numbers mirror the match 
or the mismatch of the bands of Cr and Fe at the Fermi 
level and at this particular k|| . The band structure along 
k z at k|| = (0,0.5)^ for Cr and Fe is given in Fig. §. 
For both Fe and Cr the Fermi level is aligned to be at 
zero energy. 

One can see from Fig. |(b) [k,, = (0,0.42)^] that 
\r + \ changes very rapidly around the Fermi energy. This 
means that the stationary phase approximation is not 
very well satisfied for this k|| . Therefore, the result for 
the coupling strength for X is not very reliable. For N% 
on the other hand, the variation of \r + \ and \r~\ around 
the Fermi level is smooth. The majority electrons of the 
ellipsoid are reflected strongly at the interface, because 
they have to transmit from a Cr sp-band to an Fe d- 
band (see Fig. ||). The minority electrons have smaller 
reflection amplitude \r~\ because the related Fe band 
at the Fermi level is of the same character as the Cr 
band. As soon as the energy is raised by « 0.01 Ry, 
the Fermi level minority band of Fe changes symmetry, 
and the minority electrons become completely confined 
at the Cr layer. 

Our calculations give a period of 13.77 A for N 1 . This 
period is slightly smaller than the experimental value for 
the long-period oscillation, which is somewhere between 
14 A and 18 A (see for example Refs. F Fj tnis 

same vector, the de Haas-van Alphen experiments! 2 ^ give 
a period of 15.97 A, which is within the experimental 
uncertainty. With a satisfactory period and coupling 
strength (as we will argue in the Sec. IV), the vector Ni 
is the best candidate for the origin of the long period 
oscillations of IEC. 

The very existence of the lens is due to the spin-orbit 
coupling that we have included for Cr. In the case of van- 
ishing spin-orbit parameter £ the lens touches the out- 
side surface and there is no extremal vector at this point. 
With £ = 2.84 mRy we obtained two extremal vectors L\ 
and T7 2 at k|| = (0, 0.29)^, which have periods 20.85 A 
and 17.77 A respectively. For both vectors, the factor 
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|Ai?| is not small, but in contrast to the nesting vector 
case, the curvature radius factor reg- is very small. Thus 
the overall contribution to the coupling is almost negli- 
gible. This is a somehow unwelcome result because the 
periods of L\ and L2 are very close to the experimental 
one. However, it is a result that one would anticipate 
from qualitative reasoning alone, since the phase space 
associated with the lens is small, as has been stressed by 
other authors previously!! 

The nesting vector V n corresponds to a period of about 
3.12 A. This is about 2.17 ML of Cr and agrees very well 
with previous results. For this vector, the geometrical 
weight is very big (approximately at least an order of 
magnitude bigger than most of the other Gg 's presented 
in Tabic ||) . This is as expected qualitatively, because the 
nesting feature (nearly parallel lines) leads to very large 
value for re?-. Since the reflection factor |Ai?| is not very 

small, the overall coupling strength associated with V n 
is large (37.99 mJ/m 2 ). Although this number is too big 
compared to experimental strengths, the result agrees 
with the experiment qualitatively, since it renders V n the 
best candidate for the short period oscillations (for which 
the period is measured to be 2.1 ML&Ej). Nevertheless, 
one has to acknowledge that the nesting feature requires 
a treatment that goes beyond the stationary phase ap- 
proximation used in our calculations. Specifically the 
replacement of the Fermi surface sheets by parabolas is 
questionable for this case, and one should perform an 
explicit numerical integration over ky for this particu- 
lar area of the Fermi sphere. Moreover, the coupling 
strength of the short period oscillations is very sensitive 
to roughness, as will be discussed in Sec. IV. 

B. (110) orientation 

In Table [n] we give the results of the calculation for the 
(110) orientation. Again, as in the (001) case, one can 
find, among others, extremal vectors that span the lens 
(L^, L' 2 , and L' 3 ) and vectors (for example N[) at the N 
ellipsoid. As we can see in this case, the most important 
contribution to the coupling is made by N[. Its period 
and strength are Appd to be 16.46 A (the experimental 
period is w 18 A&Ej) and 4.65 mJ/m 2 respectively. In 
Fig. [|, we give the moduli of the reflection amplitudes 
\r + 1 and \r~ | as a function of energy. The majority elec- 
trons are strongly confined in the spacer, whereas the 
minority electrons are only partially reflected at the in- 
terface. The change of |r + | and \r~\ around the Fermi 
level is slow enough, that the stationary phase approxi- 
mation is reliable. 

The lens vectors L' v L' 2 , and L' 3 have again small Kfj, 
and small overall coupling strengths. A significant con- 
tribution is made from the vector X' at k|| = (0, 0.44)^-. 
Note that we have used new x', y', z' axes so that 
2! is the growth direction. Another ellipsoid vector at 



k,| = (0.35,0.50)^ has also large J*. Anyway, these 
two JPS'g (0.89 mJ/m 2 and 1.17 mJ/m 2 ) are 5 and 4 

times smaller than that of N[. Therefore, the origin of 
the long period oscillations is attributed again to an el- 
lipsoid spanning vector. There is also a non-negligible 
strength associated with short period oscillations. This 
is for k|| = (0.0, 1.0)^, with a period of 4.70 A and a 
strength of 2.02 mJ/m 2 . To the best of our knowledge, 
no short period has been reported so far for the (110) 
orientation. 



IV. DISCUSSION OF RESULTS 

As we stressed in the previous sections, the two most 
important issues that one has to address for IEC are the 
period of the oscillation and its magnitude. For the three 
cases of (001), (110), and (211) orientation, the experi- 
mental data give a similar period. The common period 
is in favor of the argument that the large wavelength os- 
cillations are originated from a relatively isotropic part 
of the Cr Fermi surface. As we saw above, our calcu- 
lations agree with that, since the two vectors Ni and 
N[ are extremal vectors of the ellipsoid centered at N, 
which is fairly isotropic. The size of the ellipsoid, as we 
calculated it in the (110) case (period 16.46 A) or as it 
can be calculated from de Haas-van Alphen experiments 
for the (001) case (period 15.97 A), is in good agreement 
with the observed period. 

The measured coupling strengths for the (001) orien- 
tation vary from 0.6 mJ/m 2 -.ta_L6 mJ/m 2 for thick- 
nesses between 4 and 8 MLuOt* 2 !. Our own result is 
7.94 mJ/m 2 for D= 10 A, if we use Eq. (10), whereas 
if we use Eq. (3), we find that the first antiferromag- 
netic peak strength is 14.45 mJ/m 2 for D =6.7 A. These 
numbers are about 10 times the experimental ones. Pre- 
vious total energy calculations (see section 2.1 in Ref. ||) 
predicted larger amplitudes (up to 100 bigger than the 
experiment), so our result is an improvement. In a most 
recent first principle calculation, StilesEj also attributed 
the origin of the long period oscillations to N% and N[ . 
His calculated coupling strengths for these vectors are 
5.7 mJ/m 2 and 3.2 mJ/m 2 . Both numbers are in good 
agreement with ours. 

The discrepancy with the experiment can be remedied 
if the interface roughness is taken into account. The 
coupling strength in realistic samples decreases mostly 
because of the interface roughness. In the experiment, 
the Fe/Cr interfaces are not ideally flat, but they pos- 
sess steps, that lead to a random growth front of the Cr 
spacer. One can deduce from the experimental data that 
early results, obtained from rough samples grown with 
sputtering, gave smaller strengths compared to more re- 
cent ones from smooth epitaxial films. Our calculation 
can be modified in two ways in order to include the ef- 
fects of interface roughness. First, we can average the 



4 



amplitude of IEC over the variable thickness of the film. 
As is shown in Refs. |lj,[l3|, this can be done by convo- 
luting the coupling strength J(t) for a certain thickness 
t with a distribution (constant or Gaussian) over t fluc- 
tuations. Consequently, the new coupling is 

Jrou a h(D)=J2P(D,t)J(t), (11) 
t 

where D is the average thickness. Second, because of the 
roughness, the in-plane component of the wavevector is 
no longer a good quantum number. A certain ky state 
can be reflected to a range of other ky 's states around it, 
and this leads to a decrease of the coupling, or even to 
complete sweep of the oscillations, if the period is short. 
Therefore, if the sample is not prepared to be smooth, 
the 2.2 ML period associated with the nesting vector 
will not be observed. This has been verified experimen- 
tallyatl Interface roughness may also be responsible for 
the lack of the period of 4.70 A that we found in the 
(110) case. Apart from roughnesSp^other disorder effects, 
like misfit dislocations and strair£3, also can affect IEC 
in the same way. These effects however do not play any 
significant role in Fe/Cr multilayers because of the excel- 
lent match of the lattice constants (2.88 A for Cr, 2.87 A 
for Fe). 

The number that we give as the amplitude is really 
an upper bound of the actual value, since we are taking 
the absolute value of the complex number |Ai?|. This is 
demonstrated in Fig. ^, which shows vs. thickness 
D as calculated from Eqs. (3) and (10) for the span- 
ning vector N±. We use \AR\ to estimate the strength of 
IEC, because there is a certain phase arbitrariness in the 
calculation. The main source for this is the finite thick- 
ness of the Fe layers, which is not taken into account, in. 
our calculation. It has beetupredicted from theorjilja 
and verified experimentallyEil, that IEC exhibits oscilla- 
tions, not only as a function of Cr thickness, but also as 
one varies the thickness of the Fe layers. This means of 
course that the antiferromagnetic peaks occur at differ- 
ent D for different Fe widths, and so one needs to use 
the right D in each case. Moreover, for finite thickness, 
quantum well states are expected to be formed also in 
the magnetic layers and hence the coupling strength will 
be affected as well. Finally, we must point out once more 
that our approximations do not hold for very thin mul- 
tilayers, and hence the initial phase is pretty arbitrary. 

Our results indicate that the nesting vector of the el- 
lipsoid is a favorite candidate for the long period oscilla- 
tions. However, other candidates have been proposed in 
the literature as well. Early on, it was realized that the 
nesting vector can be responsible, not only for the small, 
but also for the large wavelength oscillations. Because 
of the discreteness of the thickness D, a short period 
can give rise to an effective long period as a higher har- 
monic. The phenomenon is called aliasing and one gets 
for the second harmonic a period A' = A/(A — 2) ML. 



For A « 2.2 ML the effective period is indeed large. 
However, careful measurements of the large period by 
Pierce et aln gave an effective period A' = 20.05 ML, 
much bigger than the observed period which is between 
10 to 12 ML. Moreover, the fact that the period seems to 
be independent of the orientation is not in favor of the 
aliasing scenario. Of course from the theoretical point 
of view, a conclusive analysis would be to calculate the 
coupling strength for the second harmonic and see if it is 
significant. The nesting vector though is a difficult case, 
as we mentioned in Sec. III. If the coupling strength for 
the second harmonic turns out to be comparable to the 
one of the vector Ni, then the superposition of the two 
oscillations will lead to new features in the overall IEC, 
particularly it will raise its period and also decrease its 
amplitude. 

Two other candidates for the long-period coupling are 
the lens vectors L\ and L^. In a recent paper, Li et al.c3 
have studied the confinement of electrons in the spacer 
using angle-resolved photoemission. They found sub- 
stantial confinement for states with ky in the vicinity of 
the lens. From this they concluded that the lens is the 
k-space origin of the long period oscillations. From our 
calculations we have found that the confinement at the 
lens area is big enough that one could detect quantum 
well states for this ky (the magnitude of both r + and 
r~ in the ferromagnetic configuration is w 0.5). How- 
ever, the confinement is much stronger at the ellipsoid 
as one can see in Fig. 0(a). Hence it would be useful 
to get more data for this k-space area. The confine- 
ment is one important factor for IEC, but there are also 
other factors that affect the amplitude of the oscillation. 
Most important one is the spin-asymmetry of the reflec- 
tion amplitude (as a difference of two complex numbers), 
and not individual magnitude by itself. The geometrical 
weight of the vector is a crucial factor as well, at least 
within the stationary phase approximation. For the lens 
vectors, G^-'s are very small. 

V. SUMMARY 

By using an empirical TB method (including s, p 3 , 
and d 5 orbitals) and including the spin-orbit interaction, 
we have calculated the reflection amplitudes and trans- 
mission coefficients for an Fe/Cr/Fe sandwich structure. 
Based on the spin-asymmetry of these quantities and 
the geometrical weight of the extremal vectors of the Cr 
Fermi surface, we have obtained the coupling strengths 
associated with each vector. For both (001) and (110) 
growth direction, the origin of the long period oscilla- 
tions of the exchange coupling is attributed to extremal 
vectors of the ellipsoid centered at point N. The period 
agrees satisfactorily with the observed value. The cou- 
pling strength is roughly an order of magnitude bigger 
than the experimental one, but the inclusion of inter- 
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face roughness in the calculation can improve the agree- 
ment substantially. Although the attribution of the long 
period-coupling to the aliasing effect of the nesting vec- 
tor faces several opposing arguments, one should try to 
calculate more carefully the coupling strength of the sec- 
ond harmonic of this vector. Within the approximations 
that we have used, all the other candidates for the long 
period (and in particular the lens vectors) have been 
found to have a small to negligible contribution to the 
overall coupling. 
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FIG. 1. (100) cross section of Cr Fermi surface. 



FIG. 2. Moduli of the reflection amplitudes vs. energy for 
(001) orientation, (a) k,, = (0,0.5)^ , (b) k,, = (0,0.42)^ 
(Fermi energy is aligned to be zero) . 
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FIG. 3. Band structure along line XN (ky = (0,0.5)^ for 
(a) Cr, (b) majority Fe, and (c) minority Fe. 

FIG. 4. Moduli of the reflection amplitudes vs. energy 
for (110) orientation and ky = (0.71, 0)2^ (Fermi energy is 
aligned to be zero). 

FIG. 5. Coupling strength J a vs. thickness D, as calcu- 
lated from a) Eq. (3) (solid line), b) Eq. (10) (dashed line). 



TABLE I. Results for (001) orientation. 



vector k|| [2ir/a] 


A (A) 


G a [meV] 


\AR\ 


J a (mJ/m 2 ) 




(0,0) 


5.56 


8.48 


1.01 


0.68 




(0,0) 


9.20 


160.12 


9.0 10~ 5 


2.3 10~ 3 




(0,0) 


4.18 


1.87 


0.88 


0.13 




(0,0) 


33.65 


3.40 


6.9 10" 3 


3.8 10" 3 


v n 


(0,0.25) 


3.12 


146.88 


0.40 


37.99 


u 


(0,0.29) 


20.85 


1.00 


1.7 10~ 3 


5.3 10" 4 


L 2 


(0,0.29) 


17.77 


1.39 


0.41 


0.18 


X 


(0,0.42) 


11.21 


36.99 


0.24 


2.81 


TvT 


(0.0,0.5) 


13.77 


24.06 


1.03 


7.94 


N 2 


(0.5,0.5) 


10.38 


10.40 


0.63 


1.03 




(0,1.0) 


3.42 


6.33 


0.65 


0.33 


TABLE II. Results for (110) orientation. 


vector 


ky [2tt/o] 


A (A) 


G" [meV] 


|AJ?| 


J a (mJ/m 2 ) 




(0,0) 


5.86 


18.78 


0.58 


0.86 


L[ 


(0,0.29) 


21.26 


1.37 


0.97 


0.21 


L' 2 


(0,0.29) 


17.87 


1.57 


1.11 


0.28 


L' 3 


(0.23,0) 


29.68 


3.86 


6.0 10~ 3 


3.6 10~ 3 


X' 


(0,0.44) 


9.42 


6.01 


0.93 


0.89 


W{ 


(0.71,0) 


16.46 


38.45 


1.51 


4.65 


K 


(0.35,0.50) 


12.12 


15.41 


0.24 


1.17 


Ni 


(0.71,1.0) 


10.04 


7.67 


0.25 


0.15 




(0.3,0.0) 


9.73 


12.93 


0.07 


0.30 


S 7 


(0,1.0) 


4.70 


43.31 


0.59 


2.02 
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